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ABSTRACT 

We present Gemini/GNIRS spectroscopic observations of 4 z— band [z w 7) dropout 
galaxies and VLT/XSHOOTER observations of one z— band dropout and 3 Y — band 
(z ~ 8 — 9) dropout galaxies in the Hubble Ultra Deep Field, which were selected with 
Wide Field Camera 3 imaging on the Hubble Space Telescope. We find no evidence of 
Lyman-a emission with a typical 5<7 sensitivity of 5 x 10 _18 erg cm _2 s _1 , and we use 
the upper limits on Lyman-a flux and the broad-band magnitudes to constrain the 
rest-frame equivalent widths for this line emission. Accounting for incomplete spectral 
coverage, we survey 3.0 z-band dropouts and 2.9 F-band dropouts to a Lyman-a 
rest-frame equivalent width limit > 120 A (for an unresolved emission line); for an 
equivalent width limit of 50 A the effective numbers of drop-outs surveyed fall to 1.2 
z-band drop-outs and 1.5 F-band drop-outs. A simple model where the fraction of high 
rest-frame equivalent width emitters follows the trend seen at z = 3— 6.5 is inconsistent 
with our non-detections at z = 7 — 9 at the « la level for spectrally unresolved lines, 
which may indicate that a significant neutral HI fraction in the intergalactic medium 
suppresses the Lyman-a line in z-drop and F-drop galaxies at z > 7. 

Key words: galaxies: evolution - galaxies: formation - galaxies: starburst - galaxies: 
high-redshift - ultraviolet: galaxies 



1 INTRODUCTION 

The Lyman break technique has proven to be an efficient 
tool in the selection of high-redshift candidate galaxies - 
and rest-frame UV-selected galaxies are now being regu- 
larly identified at redshifts z > 6 (e.g. Bunker et al. 2004; 
Bouwens et al. 2006; McLur<Tet al. 2010). With the ad- 
vent of the Wide Field Camera 3 (WFC3) on the Hubble 
Space Telescope (HST), the technique has been pushed to 
even higher redshifts and the number of candidates beyond 
z = 6.5 has increased to about a hundred (e.g. Bouwens et al. 
2011, McLure et al. 2011, Lorenzoni et al. 2011, Wilkins et 
al. 2011, Finkelstein et al. 2010, Wilkins et al. 2010, Bunker 
et al. 2010, Oesch et al. 2010a). Spectroscopic confirmation 
of these high-z candidate galaxies remains a very important 
task. However, the sample of convincing spectroscopically- 
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confirmed sources remains small (eg. Ono et al. 2012, Pen- 
tericci et al. 2011, Schenker et al. 2012, Vanzella et al. 2011). 
Not knowing the precise redshift of most Lyman break can- 
didates by means of spectroscopic confirmation reduces our 
confidence in inferred properties about the universe at this 
epoch, both due to redshift uncertainties on each source and 
because of the unquantifiable risk of contamination in these 
samples (e.g. see Hayes et al. 2012). For objects at such high 
redshifts, the only currently-feasible spectroscopic redshift 
diagnostic is the Lyman- a emission line, which arises from 
photoionization of HII regions by star formation. The line 
itself is resonant and sensitive to the ionization state of the 
intergalactic medium (IGM) and thus, its visibility at the 
end of cosmic reionization is expected to be reduced com- 
pared to that at lower redshifts due to the damping effect of 
an increasingly neutral IGM (Gunn & Peterson 1965; Becker 
et al. 2001). The emergence or non-emergence of Lyman-a 
may be indicative of the size of the HII regions surround- 
ing these galaxies; in order for Lyman-a emission to emerge 
from a galaxy, the ionized region surrounding the galaxy 
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needs to be sufficiently large to allow the wavelength of 
the Lyman-a photon to redshift to longer wavelengths, and 
hence become non-resonant when encountering the neutral 
IGM, thus managing to escape. 

In order to address the question of the emergence 
of Lyman-a at these high redshifts, we undertook Gem- 
ini/GNIRS and VLT/XSHOOTER spectroscopy of a sam- 
ple of 7 Lyman-break selected candidate galaxies at z > 7, 
identified as z-band and F-band dropouts with HST/WFC3 
imaging (Bunker et al. 2010, Lorenzoni et al. 2011, Wilkins 
et al. 2011) centered on and around the Hubble Ultra Deep 
Field (HUDF, Beckwith et al. 2006). We have already de- 
scribed the analysis of one of these objects (HUDF.YD3) in 
a separate paper (Bunker et al. 2012), following the reported 
detection of Lyman-a at z = 8.55 by Lehnert et al. (2010), 
so we do not discuss it again here. 

The structure of this paper is as follows. We describe 
our spectroscopic observations and data reduction in Sec- 
tion [21 and present the results of the spectroscopy in Sec- 
tion [3] We discuss our analysis in Section U and our con- 
clusions are summarized in Section [5] We adopt a ACDM 
cosmology throughout, with Qm = 0.3, Ha = 0.7 and 
Ho = 70kms _1 Mpc -1 . All magnitudes are in the AB sys- 
tem (Oke & Gunn 1983). 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 Observations with Gemini/GNIRS 

Four of these objects had initially been selected as candi- 
date z > 7 sources with HST/NICMOS imaging as they 
were undetected in the ACS images of the HUDF (Bunker et 
al. 2004) and re-confirmed with HST/WFC3 (Bunker et al. 
2010, Oesch et al. 2010a). All four of these objects have sub- 
sequently appeared in the independent selection of several 
different groups (see Wilkins et al. 2011 for a comparison). 
We also obtained spectroscopy on 5 additional candidate 
galaxies at z > 7 that were selected with i£ST/NICMOS 
but were later rejected as secure candidates by the deeper 
HST/WFC3 imaging. For this reason, we do not include 
these 5 objects in our current analysis. 

We observed the four z-band drop-out high-redshift 
galaxy candidates HUDF.zDl, HUDF.zD2, HUDF.zD3 and 
HUDF.zD4 from Bunker et al. (2010), hereafter called zDl, 
zD2, zD3 and zD4. We used the GNIRS spectrograph on the 
Gemini South telescope as part of programmes GS-2004B-Q- 
19 and GS-2005B-Q-18 (PI: A. Bunker, see also Stanway et 
al. 2004a). GNIRS is a spectrometer that can perform both 
long-slit, single-order spectroscopy in the 1.0 — 5.4/xm range 
and cross-dispersed spectroscopy in the 0.9 — 2.5/mi range. 
We used the latter cross-dispersed mode with the "high- 
dispersion" 110.5 1/mm grating and short blue camera, with 
a 7 /; -long slit. In this mode, we do not get continuous spec- 
tral coverage since not all of the spectral orders fall on the 
array. 

Three of our objects (zD2, zD3 and zD4) were observed 
using a grating angle of 21 degrees and the XD.G0525 fil- 
ter, which allowed us to include the regions 0.86-0.94/xm and 
1.005-1.095/xm in our range of covered wavelengths, corre- 
sponding to redshifts z ~ 7 for Lyman-a (as expected for 
z-band dropouts). The observations of zDl were carried out 



using two observing settings. For one setting, we used the 
same grating angle of 21 degrees but used the XD_G0507 
filter instead, whereas for the second setting we changed 
the grating angle to 23.2 degrees and used the XD.G0525 
filter. This allowed us to fill the missing wavelength range 
between 0.94/xm and 1.005/xm, as well as covering the re- 
gion 0.83/xm-0.90/xm, thus having continuous spectroscopy 
between 0.83/xm and 1.095/xm (z = 5.8 — 8.0 for Lyman-a). 

The coordinates for these four objects are taken from 
Bunker et al. (2010). In acquiring the object, a nearby star 
(with accurately determined astrometry from HST) was first 
centered in the slit and then a blind offset (less than 1') was 
executed to place the faint z-dropout galaxy in the slit. The 
Gemini telescope can execute an offset of this size to an ac- 
curacy of better than 0"£| when the same guide star is used 
throughout (as was the case for our observations). A slit 
width of r /675 was used for all observations and the data 
were acquired in a three-point dither pattern (ABC) at po- 
sitions + l / /8, 0, -1//8 along the slit long axis. We simulated 
the effect of offset errors when calculating the slit losses, and 
concluded this had at most a 10 per cent effect on the flux 
captured by the aperture. Each frame was a 900s exposure. 
In total, 9 frames were combined for each setting of zDl, 
and also for zD3 and zD4, whilst 7 frames were combined 
for zD2. The observations were conducted between 2004 and 
2006 (see Table [1]). The seeing measured from the point 
spread function (PSF) of a star in the acquisition spectrum 
was about ft' 5' Full Width at Half Maximum (FWHM). The 
seeing was monitored by the observers in real time via the 
Peripheral Wavefront Sensor star, and if the seeing degraded 
then the queue-mode observations of this programme were 
terminated and observations switched to a programme with 
a less stringent seeing requirement. In some instances during 
the observation of our target for extended periods, a reacqui- 
sition of the blind offset star was performed. We confirmed 
that the blind offset star appeared in the expected position 
(i.e. the telescope pointing had not wandered off), and we 
also checked that the seeing was again consistent (about 0"5 
on average) . Hence we are confident of the seeing used to es- 
timate the slit losses. We carried out flux calibration using 
observations of a spectrophotometric A0V telluric standard 
star which was observed on the same dates as the objects. 
From unblended spectral lines in the sky spectra we mea- 
sured a spectral resolving power of R = A/AAfwhm = 3700. 

We reduced the GNIRS spectra with IRAF using our own 
customized reduction pipeline, which is partly based on an 
existing Gemini pipeline. We first used the IRAF GNIRS task 
nvnoise to remove a vertical striping pattern that arises 
from offset bias levels. We then used nsprepare to create 
variance and data quality planes for our frames. We per- 
formed a first-pass of cosmic ray rejection using the IRAF 
tasks imcombine and ccdclip and input CCD noise and gain 
parameters. We then used the resulting rejection frames to 
mask those pixels affected by cosmic ray strikes. We used the 
same task a second time to apply a further pass of cosmic 
ray rejection, also using the data quality array to mask bad 
pixels. Flat-field generation was achieved by using the tasks 
nsr educe and nsflat on three sets of flats - each appro- 
priate for different orders: the first set being IR lamp flats, 
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the second set being short-exposure Quartz Halogen (QH) 
lamp flats and the third set being long-exposure QH lamp 
flats. Subsections of the data corresponding to each spectral 
order were then trimmed out from the original 2D spectrum 
and these were sky-subtracted to first order by subtract- 
ing the average of the other dither positions and flat-fielded 
using nsr educe. We rectified the spectra (spatially and spec- 
trally) using two-dimensional arc-spectra through a pinhole 
file, correcting for geometric distortion of the optics in the 
detector, and removed any remaining skyline residuals using 
the background task. We then added on a model of the sky- 
lines to every frame prior to combining them with imcombine 
to enable the rejection of any remaining cosmic ray strikes 
using a Poisson noise model. Finally, the sky model was 
subtracted from the resulting frame. We obtained flux cal- 
ibration using the above mentioned telluric standard star 
by means of standard IRAF techniques. The star was ex- 
tracted using the apall task. We used standard to inte- 
grate the observations of the standard stars over calibration 
bandpasses, obtain a flux using a blackbody flux distribu- 
tion model and apply an extinction correction. We then used 
the task sensf unc to determine the system sensitivity and 
finally we used the calibrate task to flux calibrate the data. 



2.2 Observations with VLT/XSHOOTER 

We observed the two F-band dropouts HUDF.YD20 (Loren- 
zoni et al. 2011) and ERS.YD2 (Lorenzoni et al. 2011.) and 
the z-band dropout P34.z.480^] (Wilkins et al. 2011) using 
the XSHOOTER spectrograph (D'Odorico et al. 2006) on 
the ESO VLT-UT2 (Kueyen) as part of programme 086.A- 
0968(B) (PI: A. Bunker). XSHOOTER is an echelle spectro- 
graph, with UV, visible and near-infrared channels obtaining 
near-continuous spectroscopy from 0.3/xm to 2.48/xm with a 
l"2-wide and ll"-long slit. HUDF.YD3 has already been 
discussed in detail in Bunker et al. (2012), so here we shall 
focus our attention on our spectroscopy of the other tar- 
gets. We dithered the observations in an ABBA sequence. 
For HUDF.YD3 and ERS.YD2, we dithered at positions +3" 
and —3" from the central coordinates along the slit long axis 
(i.e. a 'chop' size of 6"). In the case of P34.Z.4809 we dithered 
+2" 5 and — 2 / /5 from the central coordinates along the slit 
long axis (i.e. a 'chop' size of 5"). 

For ERS.YD2 we set the central coordinates to those 
from Lorenzoni et al. (2011) - see Table [2] We set the po- 
sition angle (P. A.) to be +31° (measured East of North). 
We first peaked up on a nearby star 68" 51 East and 25 / /7 
South of the desired central pointing, then did a blind 
offset. For P34.Z.4809 we set the central coordinates to 
be RA=03:33:03.765 Dec.=-27:51:20.11 (J2000) so that we 
could target both the position of P34.Z.4809 in Wilkins et al. 
(2011) - see Table[2] and UDF092y-03751196d from Bouwens 
et al. (2010) which is 0"8 away. The P.A. was set to -22.5° 
to intercept both objects. We first peaked up on a nearby 



2 This corresponds to UDFy-38135539 in Bouwens et al. (2010), 
1721y in McLure et al. (2010), z8-B115 in Yan et al. (2010), object 
125 in Finkelstein et al. (2010) 

3 This corresponds to the F-drop UDF092y-03781204 in Bouwens 
et al. 2011. 



star 17 / /98 East and 10" South of the desired central point- 
ing, and then again did a blind offset. For blind offsets of 
this size, ESO guarantee an accuracy better than 0"l where 
the guide star remains the same. We simulated the effect of 
this positional uncertainty on the fraction of the light falling 
down the slit, and concluded this was at most a 5 — 10 per 
cent effect. 

The ERS.YD2 observations were conducted in 6 ob- 
serving blocks, each including 49mins of on-source integra- 
tion and consisting of a single ABBA sequence with three 
exposures of the near-IR arm of duration 245 s at each A 
or B position. The observations were taken on the nights 
of UT 2010 December 07 and 2011 January 02, 04, 05, 11 
& 23. The vast majority of the frames were taken in good 
seeing conditions of 0" 56-0" 76 FWHM. The P34.Z.4809 ob- 
servations were conducted in 5 observing blocks, 49mins of 
which consisted of on-source integration. These were taken 
on the nights of UT 2010 October 16, 17, 19 & 28 with 
two observing blocks taken on the night of UT 2010 Octo- 
ber 17 and single observing blocks on the other nights. See- 
ing conditions were similar to those during the observations 
of ERS.YD2. All observations were taken at low airmass, 
with an average airmass of 1.16, and 83 per cent of the ob- 
servations were taken at airmasses below 1.3. At such low 
airmass, the effect of differential atmospheric dispersion is 
negligible between the red end of the optical channel and 
the near- infrared 1 — 3 fim (where we expect Lyman-a), and 
we confirmed that the alignment star was centred in the slit 
in both the optical and near-infrared spectra. Three piezo 
controlled mirrors, located in front of each arm, guarantee 
that the optical path is maintained aligned against instru- 
ment flexure and correct for differential atmospheric refrac- 
tion between the telescope guiding wavelength and each arm 
central wavelength. 

The resolving power attained for our IR-channel obser- 
vations of ERS.YD2 was R = 5000. In the case of P34.z4809 
we also considered data acquired with the optical arm as 
well as that from the near-infrared since this object is a z- 
drop, so the expected redshift for Lyman-a range extends 
down to just below 1/im. These optical data were obtained 
simultaneously with the IR data and were acquired in 5 
observing blocks, each including 49 minutes of on-source in- 
tegration. The resolving power for the optical channel was 
R = 6700. There is a small region of overlap between the 
two arms (0.994 — 1.013 /mi) and in this range we combined 
the 2D spectra using inverse-variance weighting as a func- 
tion of wavelength, which accounted for the strong variation 
of throughputs/sensitivities with wavelength in the overlap 
region and also the different readout noise characteristics of 
the two detectors. 

We used the ESO pipeline (Modigliani et al. 2010) to 
reduce our data. This pipeline applied spatial and spectral 
rectification to the spectra (which exhibited significant spa- 
tial curvature as well as a non-linear wavelength scale) by 
using the two-dimensional arc spectra through a pinhole 
mask. For the IR channel, the pipeline mapped the data 
to an output spectral scale of 1 Apix -1 and a spatial scale 
of O'.^l (from original scales of about 0.5 Apix -1 and 0"24 
respectively). For the z-drop P.34.Z.4809 we might expect 
that Lyman-a falls at the red end of the optical spectrum 
(around 0.9 /mi), hence for this object we also inspected the 
optical channel, where the data were mapped to an output 
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spectral scale of 0.4Apix _1 and a spatial scale of 0V16. In 
both channels the data were flat-fielded and cosmic rays were 
identified and masked using the algorithm of van Dokkum 
(2001). The two dither positions were subtracted to remove 
the sky to first order, and the different echelle orders were 
combined together into a continuous spectrum (taking into 
account the variation in throughput with wavelength in dif- 
ferent overlapping echelle orders) before spatially registering 
and combining the data taken at the two dither positions, 
and removing any residual sky background. Flux calibra- 
tion was achieved through observations of standard stars 
LTT3218, GD-71 and Feige 110 taken on the same nights as 
the science data. We chose the subset of the standard star 
observations that were taken in similar seeing conditions and 
at similar airmass to our science data. 



3 RESULTS 

We inspected all the 2D spectra, focusing in particular on 
the expected location for the target and the wavelength 
ranges between 0.8 and 1.2 microns where Lyman-a might 
be expected. We did this through visual inspection, includ- 
ing examining frames which had been smoothed by means 
of convolution with a Gaussian with similar FWHM to the 
spatial seeing and spectral resolution to bring up any faint 
feature. We developed a noise model, based on the Poisson 
counts of the sky background and dark current, the sensitiv- 
ity of the detector (as a function of wavelength and position 
on the array) and the readout noise of the array. Dividing 
our reduced 2D spectrum by the noise model provided a 
map of signal-to-noise ratio. We ran SExtractor (Bertin and 
Arnouts 1996) on the 2D spectrum after division by our 
noise model to identify possible emission lines which might 
have been missed through visual inspection. As a check on 
contamination by spurious sources and noise spikes, we also 
examined the negative image to confirm that there were no 
significant detections. 

We do not detect any significant line emission in any 
of the objects in our sample. To test the recoverability of 
possible line emission in our spectroscopy, we added fake 
sources of various intensity in random locations in our 2D 
spectra and checked if we would have detected these. We 
initially took an emission line with an elliptical Gaussian 
profile, with a FWHM of 200 kins" 1 in the spectral direc- 
tion, and compact spatially (to match the typical small sizes 
of the Lyman break population at z > 6) . We then convolved 
this with another Gaussian to reflect the instrumental spec- 
tral resolution and the ground based seeing (i.e. FWHM 
of about 0"5 and 2.56 A for the near-infrared spectroscopy, 
and 1.4 A for the XSHOOTER optical arm). We note that 
the effect of the strong Lyman-a forest absorption at these 
redshifts would be expected to absorb the entire blue wing 
of this line emission, so the actual line width and flux be- 
fore absorption might be twice as large. To reflect this, we 
also experimented with introducing a truncated Gaussian 
before convolution with the instrumental resolution, where 
the initial FWHM was 400kms _1 and we set the blue half 
of the profile to be zero (i.e. the FWHM is now 200kms _1 ). 
This profile produced similar recoverability statistics to the 
complete Gaussian simulations. For the XSHOOTER obser- 
vations, we find that a typical Lyman-a emission line with 



an intrinsic velocity width of about 200kms _1 would be ro- 
bustly picked up in our spectroscopy if it represented a S/N 
of 3.5-4cr (using a 5A xO'.'S aperture for the Near-IR chan- 
nel and a 4A x0'.'8 aperture for the optical channel). In the 
case of our GNIRS observations, whose data reduction is 
subject to more systematics than our XSHOOTER observa- 
tions, we impose a stricter 5cr detection threshold for a line 
with intrinsic velocity of 200kms _1 . Some of the GNIRS 
systematics arise from the low- refractive index layers of the 
ant i- reflect ion coatings on the GNIRS lenses which contain 
radioactive thorium, causing the array to be peppered with 
spikes during long exposures. Additionally, the GNAAC con- 
troller of GNIRS superimposes systematics (vertical striping 
and horizontal banding). 

In conjunction with the continuum flux inferred from 
the HST imaging, our spectroscopy is deep enough to al- 
low us to place interesting rest-frame equivalent- width (EW) 
limits on Lyman-a emission from our objects as shown in 
Figures [1] - [9] The continuum is inferred from the broad- 
band photometry from HST, from the filter above the 
Lyman-a break, i.e. the F105W F-band for the z'-drops 
at z ~ 7, and the F125W J-band for the F-drops at 
z ~ 8 — 9. We assume here a rest-frame spectral slope above 
the Lyman-a break of the form f\ oc A beta , where we adopt 
f3 = — 2 (equivalent to a spectrum flat in f u ) : as found to 
be typical of z > 6 galaxies (Stanway et al. 2005; Wilkins 
et al. 2011), including the galaxies surveyed here (Bunker et 
al. 2010). We experimented with a range of spectral slopes 
between a redder slope of f3 — —1.5 and a bluer slope of 
f3 = —2.5, and found that this only changed the inferred 
equivalent width limits by about 5 per cent (because the 
filter used to infer the continuum lies at wavelengths just 
beyond Lyman-a). At the high redshift extreme of the ex- 
pected redshift distribution (Figure ITS]), the Lyman-a break 
can encroach on the shortest wavelengths of the HST band- 
pass used to infer the continuum, and we correct for this 
assuming total absorption below the Lyman-a wavelength 
at these redshifts (the Gunn- Peterson effect). In these cases 
we also correct for possible Lyman-a line emission contami- 
nation of the broad-band magnitude, given our upper limits 
on the line flux. The significance of the upper limits on the 
flux of an emission line, given our non-detections, depends 
on the wavelength (because the sky spectrum, atmospheric 
transmission and detector sensitivity vary with wavelength) , 
and also on the spatial and spectral extent of any line emis- 
sion. 

The Lyman break galaxies at z > 6 are typically very 
compact in HST images (e.g., Oesch et al. 2010b finds r w 
0.7 kpc) and are unresolved in ground-based seeing. Hence 
we adopt a spatial extraction aperture about 1.5 times the 
seeing disk in order to maximize the signal-to-noise ratio. 
The spectral extent (i.e. velocity width) of the Lyman-a 
emission is less certain. We consider two scenarios: one where 
the line emission is unresolved at the spectral resolution of 
GNIRS or XSHOOTER (i.e. A^fwhm < 100 km s -1 ), and 
the other where the intrinsic line width is around 200 km s _1 , 
similar to that seen in some Lyman break galaxies at z ~ 
6 which have Lyman-a emission (e.g. Bunker et al. 2003, 
Stanway et al. 2004b). 

For XSHOOTER, our spectral resolution is ~ 
2.5 A FWHM in the near-infrared (the optical arm of 
XSHOOTER has higher resolution of 1.4 A), and a line with 
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Table 1. z-band dropouts targetted with Gemini/GNIRS. The RA & Dec positions are from Bunker et al. (2010). The Y- and J-band 
magnitudes quoted for zDl, zD2 and zD3 are from Wilkins et al. (2010) using the WFC3 YoQSm/^ios™ and F125W niters respectively 
whereas the Y- and J-band magnitudes quoted for zD4 are from Bunker et al. (2010). Mi600 is the absolute rest-frame UV magnitude 
around 1600A for the most probable redshift. 



Object 


R.A. (J2000) 


Dec (J2000) 


y A B-mag 


J A B-mag 


^1600 


Exp Time (hrs) 


Date 


zDl 


03:32:42.56 


-27:46:56.6 


26.71 ±0.03 


26.44 ±0.03 


-20.53 


2.25 
2.25 


27 Nov 2004 
20 Dec 2005 


zD2 


03:32:38.81 


-27:47:07.2 


27.48 ± 0.06 


26.90 ±0.04 


-20.18 


1.75 


08 Dec 2005 


zD3 


03:32:42.57 


-27:47:31.5 


27.50 ±0.07 


27.10 ±0.05 


-19.64 


2.25 


17 Dec 2005 


zD4 


03:32:39.55 


-27:47:17.5 


27.84 ±0.09 


27.34 ±0.05 


-19.79 


2.25 
2.25 


18 Dec 2005 
30 Jan 2006 



Table 2. The 3 F-band dropouts and 1 z-band dropout targeted with VLT/XSHOOTER. The RA & Dec positions are from Lorenzoni 
et al. (2011) for HUDF.YD3 and ERS.YD2, from Wilkins et al. (2011) for P34.Z.4809 and from Bouwens et al. (2011) for UDF092y- 
03751 196d. The J-band magnitudes for each object are quoted from the same respective papers using the WFC3 F125W filter. M1600 
is the absolute rest-frame UV magnitude around 1600A for the most probable redshift. 



Object 


R.A. (J2000) 


Dec (J2000) 


Mag (J) 


A^1600 


Exp Time (hrs) 


Date 




HUDF.YD3 


03:32:38.135 


-27:45:54.03 


28.18 ±0.13 


-19.12 


0.82 
1.63 
1.63 
0.82 


27 Dec 2010 

29 Dec 2010 

30 Dec 2010 

31 Dec 2010 


ERS.YD2 


03:32:02.986 


-27:43:51.95 


26.98 ±0.15 


-20.28 


0.82 
0.82 
0.82 
0.82 
0.82 
0.82 


07 Dec 2010 
02 Jan 2011 

04 Jan 2011 

05 Jan 2011 
11 Jan 2011 
23 Jan 2011 


UDF092y-03751196d 
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17 Oct 2010 
19 Oct 2010 
28 Oct 2010 


P34.Z.4809 
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26.39 ±0.03 
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intrinsic velocity width of A^fwhm = 200 km s -1 would re- 
sult in observed line widths of 6.5 — 8 A FWHM, after convo- 
lution with the spectral resolution of the instruments. Hence 
we adopt a spectral extraction width of 4 A (10 pixels) in 
the optical and 5 A (5 pixels) for the infrared channels of 
XSHOOTER, intermediate in size between the wavelength 
spread of the emission lines in our two scenarios. Spatially, 
we adopt a size of O'/S for our aperture (which is 5pix in 
the optical channel and 4pix in the infrared channel) . In the 
case of an unresolved line, we capture 95.4% of the flux with 
our aperture in our optical data and 87% of the flux in our 
Near-Infrared data. For a 200kms _1 line we capture 48.5% 
of the flux with our aperture in our optical data and 52.6% 
of the flux in our Near-Infrared data. We apply the above 
aperture corrections in computing the flux limits. 

In the case of the GNIRS observations, the spectral res- 
olution is ~ 2.56 A FWHM, and a line with an intrinsic 
velocity width of A^fwhm = 200kms _1 would correspond 
to an observed line width of « 7 A FWHM after convolution 



with the spectral resolution of the instrument. We adopt a 
spectral extent of 3.96 A for our aperture (corresponding to 5 
pixels in our pipeline output) and a spatial extraction width 
of Of! 75 (5pixels). Our square aperture (5 x 5-pixels) would 
capture 88% of the flux from a spectrally unresolved object, 
whereas for a 200kms _1 it would capture 45% of the flux. 

We checked the validity of our 2D noise model by plac- 
ing the chosen extraction apertures at random on the 2D 
spectrum normalised by the noise model, and fitting a Gaus- 
sian to a histogram of the measured fluxes within the aper- 
tures (Figure [15] shows the results of this for all spatially- 
independent apertures in one of our XSHOOTER near- 
infrared spectra). The noise distribution was well- fit by a 
Gaussian with the expected noise properties. Some excess 
power in the wings (both positive and negative) was at- 
tributable to occasional isolated hot pixels or cosmic rays 
which had survived clipping in the data reduction, or sky 
line subtraction residuals. All > 4 a events (measured in the 
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Figure 7. 4a line flux limit for our observations with the Near-IR channel of XSHOOTER. 
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Figure 8. 4a Equivalent Width limit for ERS.YD2 from observations with XSHOOTER. This plot and the next two EW upper-limit 
plots for XSHOOTER observations of HUDF.YD3 and P34.Z.4809 assume a spectrally unresolved source; a typical line with intrinsic 
velocity width of 200kms _1 would have to be about 1.5 times brighter than an unresolved line to be robustly detected in our XSHOOTER 
spectroscopy (i.e. the limits plotted here correspond to w 3a for such a spectrally resolved line). 
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Figure 9. 4a Equivalent Width limit for HUDF.YD3 from observations with XSHOOTER. 
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Figure 10. Act Equivalent Width limit for UDF092y-03751196d from observations with XSHOOTER. 
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Figure 11. Act line flux limit for our observations of P34.Z.4809 with the Optical and Near-IR channels of XSHOOTER. The region of 
the figure shown with a dotted line represents the limit obtained form data acquired with the NIR channel. 
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Figure 12. Act Equivalent Width limit for P34.Z.4809 from observations with XSHOOTER. The region of the figure shown with a dotted 
line represents the limit obtained from data acquired with the NIR channel. 
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Figure 13. This figure shows typical reduced 2D frames from all the different spectroscopic settings used for our GNIRS observations, 
together with the corresponding 2D frames showing sky emission lines. From top to bottom, (a) the 0.84-0.94/im wavelength range 
covered with the XD_G0507 filter and a grating angle of 21 (order 5), (b) the same wavelength region covered with the XD_G0525 filter 
and a grating angle of 21 (order 5), (c) the 0.95-1.03/im wavelength range covered with the G_0525 filter and a grating angle of 23.2 
(order 5), (d) the 1.005-1.095/xm wavelength range covered with the G_0507 filter and a grating angle of 21 (order 4), and (e) the same 
wavelength region covered with the XD_G0525 filter and a grating angle of 21 (order 4). Wavelength increases from bottom to top. 
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Figure 14. This figure shows typical reduced 2D frames for XSHOOTER, together with the corresponding 2D frames showing sky 
emission lines. Displayed here is the wavelength range where Lyman-a is expected to lie for z-drops and y-drops. From top to bottom, 
the first four panels show reduced data from the optical channel (here showing P34.Z.4809), and together span the 0.9/im-1.0/im wavelength 
range, each panel spanning 0.025/im with wavelength increasing from left to right. The following six panels show reduced data from 
the Near-IR channel (here showing UDF092y-03751196d) and together span the 1.0/xm-1.3/xm wavelength range, each panel spanning 
0.05/im with wavelength increasing from left to right. 
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Figure 1. Line flux and Equivalent Width limits for zDl from 
GNIRS observations. The upper panel shows the 5<j line flux limit 
probed by Order 5 of our spectroscopy using a grating angle of 
21 and the XD_G0507 filter. This filter does not have good trans- 
mission in the 0.8-1.2/xm range, which is why the sensitivity is 
not very good. The lower panel shows the 5cr Equivalent Width 
limit for the redshift range probed by the same order of the same 
spectroscopic setting. Y-mag: 26.71. This plot (and all other EW 
upper-limit plots for GNIRS) assumes a spectrally unresolved 
source; a typical line with intrinsic velocity width of 200kms — 1 
would have to be about 2 times brighter than an unresolved line 
to be robustly detected in our GNIRS spectroscopy. 



Figure 2. Line flux and Equivalent Width limits for zDl from 
GNIRS observations (continued). The upper panel shows the 5cr 
line flux limit probed by Order 5 of our spectroscopy using a 
grating angle of 23.2 and the the XD_G0525 filter. The lower 
panel shows the 5<r Equivalent Width limit for the redshift range 
probed by the same order of the same spectroscopic setting. Y- 
mag: 26.71. 



adopted apertures) were investigated and none were found 
to be consistent with emission lines. 

One of the objects targeted by our GNIRS spectroscopy, 
HUDF.zDl, has previousy been investigated by Fontana et 
al. (2010, their source G2.1408) and they observe a tenta- 
tive Lyman-a emission line at 9691. 5 A with fLya = 3.4 x 
10 -18 erg cm -2 s _1 . We looked closely at this wavelength re- 
gion in our GNIRS spectrum of this object. At this wave- 
length, our measured 1 a noise using a 3.96AxO'/75 aperture 
(5x5 pixels) is 1.01 x 10 -18 erg cm -2 s _1 . For a spectrally- 
unresolved emission line with the same flux as the Fontana 
et al. object, we would expect a 3 a detection in our GNIRS 
data. However, Fontana et al. report a marginally spectrally- 
resolved line profile of 10 A FWHM (with an instrumental 
width of 7 A FWHM for FORS2), corresponding to an in- 
trinsic velocity spread of 200kms _1 FWHM, after decon- 



volution with their instrumental resolution. At our GNIRS 
resolution we would expect such a line to have 7 A FWHM. 
For this velocity profile and line flux, we would only expect 
a 1.7 a signal within our 5x5 pixel aperture. Hence our non- 
detection does not rule out the Fontana et al. line detection, 
particularly if it has significant velocity extent. We note that 
although HUDF.zDl is well resolved in the WFC3 imaging 
and indeed comprises two distinct components separated by 
2 kpc (0"37) and with half-light radii of 0.5 - 0.8 kpc (Oesch 
et al. 2010b), in ground-based seeing this galaxy pair is un- 
resolved. We have included the effects of the finite source 
size (before seeing) in our calculation of the slit losses. We 
are currently analysing recently obtained VLT/FORS2 data, 
which will shed further light on this interesting target. 
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Figure 3. Line flux and Equivalent Width limits for zDl from 
GNIRS observations (continued). The upper panel shows the 5cr 
line flux limit probed by Order 4 of our spectroscopy using a 
grating angle of 21 and the XD_G0507 filter. The lower panel 
shows the 5<r Equivalent Width limit for the redshift range probed 
by the same order of the same spectroscopic setting. Y-mag: 26.71. 



4 ANALYSIS 

4.1 The EW distribution and Constraints on the 
Neutral Fraction Xhi 

Given the fact that we detect no Lyman-a emission in any of 
our spectra, what can we say about the strength of Lyman- 
a emission in z > 7 galaxies? The broad-band filters used 
to colour-select candidates with the Lyman break technique 
introduce a selection function on the redshift. The expected 
redshift distribution used for our sample is taken from the 
simulations described in Wilkins et al. (2011) at z = 7 for 
the z-band dropouts and Lorenzoni et al. (2011) at z = 8 for 
the y~band dropouts, which are reproduced in Fig |18l We 
compute the probability of recovering a galaxy as a function 
of redshift for different rest-frame UV absolute magnitudes 
around A res t = 1600 A (Mi6oo)- For each spectroscopic tar- 
get, we calculated Mi6oo so that for every object we could 
determine the most appropriate scenario to use from the 
simulations in Fig |18l choosing the relevant curve for a given 
M1600 and field (eg. HUDF, P34 etc). These curves give us 
the probability of recovering a galaxy at a given redshift. 
We then considered different thresholds on the rest- 



Figure 4. 5cr Equivalent Width limit for zD2 from GNIRS ob- 
servations. The upper panel shows the redshift range probed by 
Order 5 of our spectroscopy using a grating angle of 21 and the 
XD_G0525 filter. The lower panel shows the redshift range probed 
by Order 4. Y-mag: 27.48. 



frame EW (50 A, 80 A and 120 A) for each of our targets. 
For each particular target and chosen EW threshold, we 
computed the fraction of our spectroscopy that had EW 
limits lower than our chosen threshold (i.e. the fraction of 
the spectrum where EW U p P eriimit < EWthreshoid, which we 
call FracEw<thres, weighted by the redshift probability dis- 
tribution for the dropout galaxy). 

By considering the likelihood function for a galaxy to 
be lying at a particular redshift, we computed the fractional 
probability, Frac 2 , that a galaxy drawn from the dropout 
sample would fall within the spectral coverage of that spec- 
trograph setup (tabulated in column 3 of Tabled). Multiply- 
ing Frac 2 by FracEw<thres (found in columns 5-7 of Table [3]) 
and computing the sum over all the galaxies observed gives 
us the effective total number of galaxies, iV e ff, where our 
sensitivity is greater than our chosen EW threshold. These 
are tabulated in Table [5] 

Given that we do not detect any Lyman- a emission in 
our data we now consider what scenarios of EW evolution 
from lower redshift we can rule out. From Poisson statistics, 
if a given model predicts that on average A ex galaxies are 
expected to be detected in the survey, the probability f n of 
detecting n galaxies is given by 
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Figure 5. 5cr Equivalent Width limit for zD3 from GNIRS ob- 
servations. The upper panel shows the redshift range probed by 
Order 5 of our spectroscopy using a grating angle of 21 and the 
XD_G0525 filter. The lower panel shows the redshift range probed 
by Order 4. Y-mag: 27.5. 



Figure 6. 5cr Equivalent Width limit for zD4 from GNIRS ob- 
servations. The upper panel shows the redshift range probed by 
Order 5 of our spectroscopy using a grating angle of 21 and the 
XD_G0525 filter. The lower panel shows the redshift range probed 
by Order 4. Y-mag: 27.84. 



_ (A ex ) n e- A - 
/n " n\ 

We detect no galaxies in our survey (n = 0), and are able 
to reject models which predict A ex detections at the (1 — 
f n ) level; i.e. a model which predicts 1 galaxy is rejected 
at the 63 per cent level (roughly corresponding to 1 a for a 
Normal distribution) , and a model which predicts 3 galaxies 
are detected is rejected with 95% confidence (corresponding 
to about 2a for a Normal distribution). 

Specifically, we can rule out at the 63 percent (~ 1 a) 
level any scenario that predicts a fraction of galaxies with 
Lyman-a emission above the threshold EW, of X^ya > 
1/^Veff (taking A ex = XLya x N e s = 1 in the above equa- 
tion) . Here we compare our results with the work of Stark et 
al. (2010), building on the previous work of Stanway et al. 
(2007) at z = 6 and Shapley et al. (2003) at z = 3. Stark et 
al. (2010) determined the fraction of dropout galaxies with 
EW > 75 A at lower redshift (z = 4 — 6.5), and showed that 
over this redshift range the fraction of strong Lyman-a emit- 
ters increased with increasing redshift. The four data points 
with error bars in Figure [T71 are from Stark et al. (2010) and 
show the fraction of strong Lyman-a emitters at different 
redshifts in their spectroscopic dropout sample, which cov- 



ers a range in UV luminosities of —19.5 > Muv > —20.5, 
similar to the range of luminosities of our higher-redshift 
sample presented here. We also add a point at z = 3 from 
Shapley et al. (2003), who found that 29 of 957 [/-band 
drop-out galaxies had EW > 80 A. We do a simple linear 
extrapolation in redshift of this trend out to z = 8.5 (dotted 
line), where we are sensitive with the y-drops. This sce- 
nario would correspond to constant evolution (linear with 
redshift) in the intrinsic Lyman-a EW distribution coupled 
with no evolution in the neutral fraction from lower red- 
shifts, with xhi — 0. The three upper-limit arrows in Figure 
[I7l show the constraints derived from our observations; the 
tail and head of the arrows represent the expected fraction of 
Lyman-a emitting galaxies with EW< 75 A and EW< 120 A 
respectively. The arrow at z = 7 shows the constraint we get 
from our z-drops, whereas the one at z = 8.5 shows the con- 
straint from our y-drops. The mean redshift for 2;-drops and 
y-drops was derived from the simulations by Wilkins et al. 
(2011) and Lorenzoni et al. (2011). The arrow in the middle 
at z = 7.8 is the constraint obtained by considering the z- 
drops and y-drops together, at a mean redshift of z = 7.8. 
For a spectrally-unresolved line, we would expect on aver- 
age 1.2 galaxies to be detected in our combined sample if the 
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Figure 15. This plot for one of our XSHOOTER Near-InfraRed 
spectra shows the histogram of measured counts in independent 
5x4 pixel apertures. The core distribution is well- fit by a Gaussian 
with the expected noise properties. The non-Gaussian extremes 
of the distribution are due to a small number of data reduction 
artefacts. 

extrapolated evolution of the Lyman-a fraction holds (i.e. 
Xl V u ~ 0.3 at z = 7.8); we do not detect any galaxies, and 
hence the hypothesis is mildly inconsistent with our upper 
limits at the 70 per cent level. If the emission lines have 
intrinsic velocity widths of ~ 200 km s -1 , then we would ex- 
pect only 0.64 galaxies to have detectable Lyman-a emission 
in our sample adopting the extrapolated XLya, so we cannot 
rule out this scenario with our current data (it is formally 
inconsistent with our upper limit only at the 47 per cent 
level) . 

4.2 Brief review of and comparison with other 
studies 

Efforts to obtain spectroscopic detection of Lyman-a emis- 
sion at z ~ 7 have also been made by other groups. In this 
section, we summarize the results of these other studies and 
compare and contrast our own results. 

Fontana et al. (2010) used VLT/FORS2 to observe 7 
Lyman-break galaxy candidates selected in the GOODS-S 
field from Hawk-I/VLT and WFC 3 / HST imaging, and other 
than the one tentative Lyman-a emission line in G2_1408 
(HUDF.zDl in the catalogue of Bunker et al. 2010) which 
we also target and has been discussed above, they detect no 
other Lyman-a emission in the rest of the sample. 

As part of the same survey, Vanzella et al. (2011) re- 
port the detection of Lyman-a emission in 2 objects in the 
BDF4 field (Lehnert & Bremer 2003). The results of the fi- 
nal sample were discussed in Pentericci et al. (2011), who 
observed 20 z-drop galaxies with VLT/FORS2 and detected 
5 Lyman-a lines in their sample. Adopting simulation tech- 
niques discussed in Fontana et al. (2010) they report that, 
on the basis of observations made at lower redshift, the prob- 
ability of detecting only 5 galaxies in their sample is below 
2%, indicating a declining fraction of strong Lyman-a emit- 
ters at higher redshift s, consistent with our results. 

Schenker et al. (2012) presented Keck LRIS observa- 




Figure 16. The calibrated GNIRS spectrum, with the location of 
HUDF.zDl and the expected wavelength of the tentative Lyman- 
ol emission reported by Fontana et al. (2010) marked with a white 
circle. Wavelength increases from left to right, and we show the 
95 A either side of 9691.5 A. From top to bottom: (a) the re- 
duced data. The three vertical lines of higher noise are due to 
night sky emission lines; (b) the reduced data convolved with 
an elliptical Gaussian with a of 1.4 pixels spatially and 2.2 pixels 
spectrally, matching the profile of a Gaussian emission line with 
FWHM of 0"5 and 2.6 A. (c) a fake source with the same line 
flux (3.4 x 10 -18 erg cm -2 s _1 ) and wavelength as the Fontana 
et al. (2010) line added into the frame. The resulting frame has 
been Gaussian smoothed. We assume a spatially and spectrally 
unresolved source. If the line is spectrally resolved, the S/N would 
be lower - panel (d) shows the expected 2D Gaussian-smoothed 
spectrum for an emission line with an intrinsic velocity width 
of 200kms _1 (FWHM) and the same line flux as reported in 
Fontana et al. (2010). 

tions of 19 sources with photometric redshifts lying in the 
range 6.3 < z < 8.8. They reported two convincing Lyman- 
a detections (ERS 8496 at z = 6.441 and A1703_zD6 at 
z = 7.045), and a marginal detection at z = 6.905 for 
HUDF09_1596. For their discussion, they added the 7 ob- 
jects discussed by Fontana et al. (2010) such that they could 
carry out an analysis over a larger sample of objects. They 
conclude that 24 of their targets (out of 26) have spectral 
coverage over more than half the likely wavelength range for 
Lyman-a (given the photometric redshifts using the EAZY 
code) . From simulations they show that out of the combined 
26 targets observed by them and Fontana et al. (2010), they 
should have detected 7-8 emission lines rather than just 2, 
again arguing for a steep decline in the Lyman-a fraction at 
z > 6.3. 
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Figure 17. Our upper limits on the fraction of high rest-frame 
equivalent Lyman-a emission at z > 7 are shown for the z-drops 
[z = 7), the y-drops (z = 8.5) and our complete sample (with 
mean z = 7.8). The head of each arrow is the 120 A limit and 
the tail of each arrow is the 75 A limit. The diamond symbols 
are results obtained at lower redshift by Shapley et al. (2003) at 
z = 3 and Stark et al. (2010) at z = 4 — 6.5. For comparison, 
we extrapolate the low-redshift trend to higher redshifts (dotted 
line). Our upper limits appear inconsistent with this extrapolation 
at the 1 a level, perhaps indicating that the IGM neutral fraction 
Xhi > at z > 7. The upper figure shows our constraints when 
considering an unresolved line, whereas the lower figure considers 
a line with an intrinsic velocity width of 200km/s. 



Ono et al. (2012) carried out observations of 11 z-drops 
in the SDF and GOODS-N using Keck/DEIMOS, detecting 
Lyman-a emission for 3 objects in their sample, one of which 
had already been detected at a lower level of significance by 
lye et al. (2006). 

All of the objects discussed above were brighter in the 
y-band than any of our spectroscopic targets. We checked 
whether we would have detected Lyman-a emission from 
any of these objects in our own spectroscopy if we had tar- 
geted them. To perform these comparisons we first consid- 
ered an appropriate line profile for each individual object; 
in each case, the emission was assumed to be spatially un- 
resolved, whereas for the spectral extent of the lines we 
adopted the published values and convolved these with the 
spectral resolution of our instrument^]. We chose our aper- 
ture to be about twice the size of the adopted emission 



4 A1703_zD6 was unresolved in the observations of Schenker et 
al. using NIRSPEC (AAfwhm = 6.5A), so here we are assuming 
that this object would be unresolved in our observations with 
GNIRS and XSHOOTER. 




Figure 18. The probability of recovering a galaxy in the sim- 
ulations described in Wilkins et al. (2011) and Lorenzoni et al. 
(2011) for different fields as a function of redshift for several dif- 
ferent absolute rest-frame Mi6oo magnitudes. In this paper, we 
use the results of these simulations to provide the expected red- 
shift distributions. The upper two figures show the simulations 
for z-band dropouts whereas the lower two figures are for y-band 
dropouts. In each case, the mean redshift is denoted by a dot. 



line profile (spatially and spectrally) and performed aper- 
ture corrections to allow for any loss in flux in our aper- 
tures. We report the results of these comparisons in Table 
[3] In summary, we conclude that nearly all objects with the 
line fluxes detected by other groups would have been clearly 
detected in our spectroscopy^. However, we note that the 
detections presented by other groups are typically for ob- 
jects with brighter broad-band magnitudes, and hence the 
equivalent widths are small. Only 4 of the 9 objects in Ta- 
ble [3] have rest-frame equivalent widths > 50 A, where our 
observations are sensitive over an appreciable redshift range. 
The others have smaller rest-frame equivalent widths, which 
would be probed only by our spectra of P34.Z.4809 and zD2 
(around 1.05/xm). This suggests that whilst there might in- 
deed be no Lyman-a emission escaping the objects which 
we targetted in our spectroscopy, the negative results could 
also be due the the possibility of the emission being fainter 
than our detection limits, particularly given that our tar- 
getted objects are fainter than any of the other objects in 
literature discussed above. It has been suggested, however, 
that the fraction of high equivalent- width Lyman-a emitters 
is larger for faint galaxies (Stark et al. 2010, 2011), which 
means that our derived upper limits on the equivalent width 



5 As can be seen from Table [3] the only two objects which we 
would not have detected in either GNIRS or XSHOOTER are 
NTTDF-6345 and NTTDF-474, whose Lyman-a flux is mea- 
sured to be quite low. ERS 8496 (which also has a relatively 
low measured Lyman-a flux) would have been undetected in our 
GNIRS observations, but would have been clearly detected in 
our XSHOOTER data. BDF-3299, on the other hand, would 
have been detected in our GNIRS observations, but not in our 
XSHOOTER data. 
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Table 3. This table shows the significance levels at which objects discussed in literature would have been detected in our spectroscopy 
if we had targetted them with Gemini/GNIRS or VLT/XSHOOTER. 



Object 


Redshift 


Mag (Y AB ) 


Lyman-a flux 
(erg cm" 2 s _1 ) 


EW / A 


GNIRS a 


XSHOOTER a 


ERS 8496 


6.441 ±0.002 


N/A 


0.91 ± 1.4 x 10- 17 


69 ±10 


1.52 


6.5 


A1703_zD6 


7.045 


N/A 


2.84 ±5.3 x 10- 17 


65 ± 12 


18.93 


16.7 


BDF-521 


7.008 ±0.002 


25.86 


1.6 ±0.16 x 10~ 17 


64±£° 


6.67 


5.93 


BDF-3299 


7.109 ±0.002 


26.15 


1.2 ±0.14 x 10- 17 


501- 


4.80 


2.00 


GN-108036 


7.213 


25.50 


2.5 x 10- 17 


33 


5.00 


3.57 


SDF-63544 


6.965 


25.10 


2.8 x 10- 17 


43 


8.62 


5.19 


SDF-46975 


6.844 


25.20 


2.7 x 10~ 17 


43 


7.20 


13.50 


NTTDF-6345 


6.701 


25.46 


0.72 x 10- 17 


15 ±3 


0.48 


2.4 


NTTDF-474 


6.623 


26.50 


0.32 x 10~ 17 


16 ±5 


0.64 


1.78 



of these objects, shown in figures [Tl - [T2l may offer stronger 
constraints on Xhi- 



5 CONCLUSIONS 

We have presented spectroscopic observations with GEM- 
INI/GNIRS and VLT/XSHOOTER of a sample of z > 7 
candidate galaxies and we fail to detect significant Lyman- 
a emission from any of them. This is consistent with the 
fraction of high rest-frame equivalent width Lyman-a emit- 
ters dropping at z > 7, as would be expected if the neutral 
HI fraction was greater at these epochs. We have also inves- 
tigated a tentative emission line published by Font ana et al. 
(2010) in HUDF.zDl (from the catalogue of Bunker et al. 
2010) and our analysis does not confirm the presence of this 
line, although we do not rule out the possibility of it being 
real, especially if it has considerable velocity extent. 

Given the lack of Lyman-a emission in our spectroscopy 
in conjunction with the continuum flux derived from HST 
imaging of these objects, we derived upper limits on the 
rest-frame equivalent width of our objects. Extrapolating 
the Lyman-a fraction observed at lower redshifts by Stark 
et al. (2010) and Shapley et al. (2003), our lack of Lyman-a 
detection rules out at a level of 1 a (70%), for spectrally un- 
resolved lines, the scenario in which the Lyman-a fraction 
evolves with the same trend found at lower redshifts. The 
limits are weaker if the lines have significant velocity width 
extent. A diminished Lyman-a fraction at higher redshift 
is consistent with other published studies. This attenuation 
in the Lyman-a fraction can be attributed either to physi- 
cal evolution of the galaxies or, more likely, an increase in 
the neutral fraction of hydrogen at z > 7, i.e. these obser- 
vations can most likely be interpreted as implying that the 
neutral fraction at z ~ 8 can be ruled out as being xhi = 
at a level of 1 a. Larger- number statistics are required to 
confirm this hypothesis at a higher level of significance. To 
this end, we have undertaken spectroscopy on a large sam- 
ple of z-band and F-band dropouts with VLT/FORS2 and 
SUBARU/MOIRCS (Caruana et al. in prep). 
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Table 4. This table for all objects targeted by our spectroscopy shows the redshift range spanned by our data for Lyman-Q! in Column 
2. Column 3 shows the fractional probability that a galaxy drawn from the dropout sample would fall within the spectral coverage of 
that particular spectrograph setup. Column 4 gives the median EW for each object for the most probable redshift range (see Figure [TH)) . 
The remaining three columns tabulate the fraction of our spectroscopy that has EW limits lower than our chosen threshold (50 A, 75 A 
and 120 A respectively.) The figures in this table are for an unresolved line. 



Object 


z- Range 


spanned 


Frac 2 


Median EW 


FraC EW<50A 


FraC EW<75A 


FraC EW<i20A 




by data (for Ly-a) 












zDl 


5.82 ■ 


- 6.40 


0.01 


54.69 





0.4583 


1.0000 




6.08 


- 6.74 


0.24 


243.2 













6.79 


- 7.45 


0.50 


29.35 


0.8349 


0.9415 


0.9990 




7.25 


- 8.05 


0.36 


85.29 


Q 


u.oozu 


u .u zyy 


zD2 


6.09 


- 6.73 


0.24 


113.5 


Q 


Q 


0.7697 




7.25 


- 8.00 


0.33 


76.10 


0.285052 


0.4109 


0.6883 


zD3 


6.09 


- 6.73 


0.24 


120.4 








0.7262 




7.25 


- 8.00 


0.33 


80.72 





0.3601 


0.6413 


zD4 


6.09 


- 6.73 


0.24 


150.2 








0.1869 




7.25 


- 8.00 


0.33 


100.8 





0.0689 


0.5326 
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3.60 - 7.40 (Optical) 


0.70 


19.45 
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0.9188 


0.9856 
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0.8191 


1.2917 


1.7287 


ERS.YD2 
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0.6529 


1.1551 


1.4293 


HUDF.YD3 


7.42 - 19.40 (NIR) 


0.98 


241.6 


0.0550 


0.1881 


0.6516 


UDF092y-03751196d 


7.42 - 19.40 (NIR) 


0.996 


28.64 


0.8400 


0.9374 


0.9872 



Table 5. This table shows the total effective number of sampled galaxies with an EW upper limit lower than a set threshold. We present 
these figures separately for z-drops, Y-drops and both z-drops and Y-drops combined. The figures in this table are for an unresolved 
line. 





Average Redshift 




A^eff = Yl Frac * X FraC E W<thres 








EW thres = 50A 


EW thres = 75A 


EW thres = 120A 


z-drops 


7.0 


1.242 


1.903 


2.963 


Y-drops 


8.5 


1.452 


2.111 


2.851 


z-drops & Y-drops 
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2.694 


4.014 


5.814 
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Table 6. This table for all objects targeted by our spectroscopy shows the redshift range spanned by our data for Lyman-a in Column 
2. Column 3 shows the fractional probability that a galaxy drawn from the dropout sample would fall within the spectral coverage of 
that particular spectrograph setup. Column 4 gives the median EW for each object for the most probable redshift range (see Figure [TH)) . 
The remaining three columns tabulate the fraction of our spectroscopy that has EW limits lower than our chosen threshold (50 A, 75 A 
and 120 A respectively.) The figures in this table are for a 200km/s line. 



Object z- Range spanned 

by data (for Ly-a) 



Frac 2 



Median EW Frac 



Frac 



Frac 



zDl 


5.82 - 


6.40 


0.01 


107.4 













6.08 - 


6.74 


0.24 


476.2 













6.79 - 


7.45 


0.50 


69.38 


0.1958 


0.5460 


0.8672 




7.25 - 


8.05 


0.36 


167.3 








0.1016 


zD2 


6.09 - 


6.73 


0.24 


222.5 













7.25 - 


8.00 


0.33 


149.3 








0.2177 


zD3 


6.09 - 


6.73 


0.24 


235.9 













7.25 - 


8.00 


0.33 


158.4 





01 


0.1062 


zD4 


6.09 - 


6.73 


0.24 


294.4 













7.25 - 


8.00 


0.33 


197.7 











P34.Z.4809 


3.60 - 7.40 (Optical) 


0.70 


38.87 


0.5759 


0.6811 


0.8617 




7.42 - 19.40 (NIR) 


0.30 


58.65 


0.2650 


0.7361 


1.2627 


ERS.YD2 


7.42 - 19.40 (NIR) 


0.86 


112.4 


0.1583 


0.4344 


1.0600 


HUDF.YD3 


7.42 - 19.40 (NIR) 


0.98 


400.3 








0.0431 


UDF092y-03751196d 


7.42 - 19.40 (NIR) 


0.996 


50.63 


0.5331 


0.7838 


0.9181 



Table 7. This table shows the total effective number of sampled galaxies with an EW upper limit lower than a set threshold. We present 
these figures separately for z-drops, Y-drops and both z-drops and Y-drops combined. The figures in this table are for a 200km/s line. 





Average Redshift 




A^eff = Yl Frac * X FraC E W<thres 








EW thres = 50A 


EW thres = 75A 


EW thres = 120 A 


z-drops 


7.0 


0.581 


0.971 


1.559 


Y-drops 


8.5 


0.667 


1.154 


1.868 


z-drops & Y-drops 


7.8 


1.248 


2.125 


3.427 
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